Some Results on Cylindrical Cavity Resonators 

By J. P. KINZER and I. G. WILSON 

Certain hitherto unpublished theoretical results on cylindrical cavity reson- 
ators are derived. These are: an approximation formula for the total number 
of resonances in a circular cylinder; conditions to yield the minimum volume cir- 
cular cylinder for an assigned Q; limitation of the frequency range of a tunable 
circular cy Under as set by ambiguity; resonant frequencies of the elliptic cylinder; 
resonant frequencies and Q of a coaxial resonator in its higher modes; and a brief 
discussion of fins in a circular cylinder. 

The essential results are condensed in a number of new tables and graphs. 

Introduction 

THE subject of wave guides and the closely allied cavity resonators was 
of considerable interest even prior to 1942, as shown in the bibliography. 
It is believed that this bibliography includes virtually everything published 
up to the end of 1942. During the war, many applications of cavity reso- 
nators were made. Among these was the use of a tunable circular cy Under 
cavity in the TE 01m mode as a radar test set; this has been treated in pre- 
vious papers. 1,2 During this development, a number of new theoretical 
results were obtained; some of these have been published. 2 Here we give 
the derivation of these results together with a number of others not previ- 
ously disclosed. 

In the interests of brevity, an effort has been made to eliminate all 
material already published. For this reason, the topics are rather discon- 
nected, and it is also assumed that the reader has an adequate background 
in the subject, such as may be obtained from a study of references 3 to 7 
of the bibliography, or a text such as Sarbacher and Edson. 8 

A convenient reference and starting point is afforded by Fig. 1, taken from 
the Wilson, Schramm, Kinzer paper. 2 This figure also explains most 
of the notation used herein. 
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1. Approximation formula for number of resonances in a circular cylin- 
drical cavity resonator. 

2. Conditions for minimum volume for an assigned Q. 
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SOME RESULTS ON CYLINDRICAL CAVITY RESONATORS 411 

3. Limitation of frequency range of a tunable cavity in the TE Oln mode 
as set by ambiguity. 

4. Resonant frequencies of an elliptic cylinder. 

5. Resonant frequencies and Q of higher order modes of a coaxial reso- 
nator. 

6. Fins in a circular cylinder. 

Approximation Formula for Number of Resonances in a 
Circular Cylinder 

From Fig. 1, the resonant frequencies of the cylindrical cavity are ob- 
tained from the equation: 

«*-©'+ (5)' (1 > 

*n which r is written in place of r/ m , to simplify the equations. The dis- 
tribution of the resonant frequencies, starting with the lowest, can be 
approximated by a continuous function 

N = F(f ) m G(X ) 

where N represents the total number of resonances up to a frequency /o 
or a wavelength Xo . This is bound to be an approximation, since the true 
function F is discontinuous (or stepped) by virtue of the resonances being a 
series of discrete values. For practical purposes, if F fits the stepped curve 
so that the steps fluctuate above and below F, it will be a useful approxi- 
mation. 

Derivation of such a formula as applied to the acoustic resonances of a 
rectangular box has recently been a subject of investigation by Bolt 9 and 
Maa. 10 Only slight modifications of their method need be made to apply 
to the present situation. 



Multiply (1) thru by 



(?)• 
(*)' 



I iran' 2 

= r +{ u 



Hence, if a point ( r, — J is plotted on the XY plane the distance from the 

origin to this point will be — and hence a measure of the resonant fre- 

c 

quency. If all such points are plotted, they will form a lattice represent- 
ing all the possible modes of resonance. The problem, then, is to find the 

number of lattice points in a quadrant of a circle with radius, R = ^-° . 
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The values of the Bessel zero, r, are not evenly spaced along the X axis; 
indeed the density, or number per unit distance, increases as r increases. 
Let the density be p(x). Then the problem becomes one of finding the 
weight of a quadrant of material whose density varies as p(x). 

Suppose the expression for M, the number of zeros r, less than some value 
x, is of the form 

M = Ax 2 + Bx 
whence, by differentiation, 

p(x) = 2 Ax + B. 

The weight, W, of the quadrant of a circle of radius R is then, by integra- 
tion, 

W = \ AR Z + - t BR 2 . 
3 4 

IT 2LW 

Since there are — lattice points per unit distance along the Y axis, 

ira na 

is apparently the total number of points in the quadrant. However, there 

are two small corrections to consider. First is that in this procedure a 

lattice point is represented by an area and for the points along the X axis 

7Tfl 

half the area, i.e., a strip — wide lying in the adjacent quadrant, has been 

omitted. Second is that the restriction n > for TE modes eliminates 
half the points along the X axis. As it happens, these corrections just 
cancel each other. Thus we have 



in which 



16wA V ,tBjS 

N = T xS + 2 xl 



4 Jo 



From a tabulation 11 of the first 180 values of r, the empirical values A = 
0.262, B = were obtained. This gives 

V 
N = 4.39 ri . 

Ao 

Subsequently, from an analysis of over a thousand modes in a "square 
cylinder" (a = L), Dr. Alfredo Bafios, formerly of M.I.T. Radiation Lab- 
oratory, has calculated the empirical formula 

N = 4.38 -3 + 0.089 3 (2) 

Ao Ao 
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from which A = 0.262, B = 0.057. These values give better agreement 

with the 180 tabulated values of r. 

There is a two-fold degeneracy in a circular cylinder for modes with 

/ > 0, which is removed, for example, when the cylinder is made elliptical. 

The total number of modes, then, counting degeneracies twice, is about 2N, 

which brings (2) in line with the general result that, in any cavity resonator, 

St V 
the total number of modes is of the order — 75 . 

3 Ao 

Minimum Volume of Circular Cylinder for Assigned Q 

In practical applications of resonant cavities, the conditions of operation 
may require high values of Q which can be attained only by the use of high 
order modes. The total number of modes, most of which are undesired, 
can then be reduced only by making the cavity volume as small as possible, 
consistent with meeting the requirement on Q. 

It will be shown that, for a cylinder, operation in the TE Olw mode very 
probably gives the smallest volume for an assigned Q. 

Statement of Problem 

When the relative proportions (the shape) of a cavity and the mode of 
oscillation are fixed, both the Q and the volume, V, of the cavity are func- 
tions of the operating wavelength, X. Since we are primarily interested 
in the relationship between Q and V, with X fixed, some simplification can 
be made by eliminating X as a parameter. This may be done by a change of 

8 V 

variables to Q - and — , respectively; to simplify the typography, these 

A A 

quantities will be denoted by single symbols: 

We are, consequently, interested in the following specific problem: 
In a circular cylindrical resonator, which is the optimum mode 

family and what is the corresponding shape to obtain the smallest 

value of W for a preassigned value of P? 

A rigorous solution cannot be obtained by the methods of elementary 
calculus, since P is not a continuous function of the mode of oscillation. 
However, a possible procedure is to assume continuity, and examine the 
relation between P and W under this assumption. If sufficiently positive 
results are obtained, the conclusions may then be carried over to the dis- 
continuous (i.e., the physical) case with reasonable assurance that, except 
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perhaps for special values, the correct answer is obtained. We proceed on 
this basis. 

Solution 

To permit a more coherent presentation of the arguments, only their 
general outline follows. More mathematical details are given later. 

We start with the formulas for Q - (= P) as given in Fig. 1. 

A 

The first operation is to show that, under comparable conditions, i.e., 
X, r, n fixed, the TE Omn modes give the highest values of P. That this is 
plausible can be seen in a general manner from the equations as they stand. 
For the TE modes, if I = 0, the numerator of the fraction is largest. Also, 
P simplifies, and the denominator roughly reduces the expression in square 
brackets to the 1/2 power. Now compare this expression with those for 
the TM modes. That for the TM modes (w > 0) is smaller because of the 
factor (1 + R) in the denominator. Finally, that for the TM modes (w = 
0) is still smaller, because 1 < (1 + p^R 2 ) 1 '-. 

This leaves only the TE Omn modes to be considered, and the next step 
is to show that m = 1 is the most favorable value. Since the relation be- 
tween P and W is cdmplicated, a parameter <p is introduced, with <p defined 
by 

tan tp = pR. (3) 

The resulting parametric equations are: 

P = - n : (4) 

lir 3 .1.3 

cos <p + - sm ip 

P 

47r cos <p sin <p 

For each of the discrete values of r and n (n is related to p) then, plots 
of P vs W can be prepared as shown in Fig. 2 for the TE Olw modes. 

Inspection of Fig. 2 shows that the best value of Q does not correspond 
to a minimum of W or a maximum of P for a given value of n, but rather to 
a point on the "envelope" of the curves. To get the envelope, we assume 
p to be continuous and proceed in the standard manner. It turns out that, 
by solving (4) for p in terms of P, r and tp, substituting the resulting ex- 

pression in W, and setting — = an equation is obtained which, when 

solved for <p, gives the values of <p which lie on the envelope. 
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dW 
We next substitute this expression for ip in W and calculate -— assuming 

now that r is continuous, and find that W has no minimum. Practically ,_ 

this means that the smallest value of r should be used, i.e., the TE 01» mode. _ . 

Finally, since from Fig. 2 it is seen that the envelope is reasonably smooth 

for values of Q ~ > 1, the expression for <p derived on the assumption of 

A -• 

continuous p is used to obtain a simple relation of great utility in practical 
cavity design. 

Details of solutioti 

In (3), since R must be finite for a physical cylinder, < tan <p < oo , 
< sin (p < 1, and < cos <p < 1. Hence we may always divide by 
sin ip or cos <p. Note that <p ranges between 0° and 90°. 

From Fig. 1, 

2r(l + p 2 R 2 ) m 



k = 



whence 



k sin (p = 



= 2 Jr3 (6) 



a 



k cos <p = — . (7) 



a 



We define W by: 



W m Y. = r ± i- . (8) 

w X 3 4/2 87T 3 



Substituting (6) and (7) in (8), 

w = S -tAjt- • (50 

4?r cos v 5 sin <P 

Substitution of (3) into the expression for Q- (= P) for the TE modes as 
given in Fig. 1 yields, after some manipulation 



f_ 1 - Wr) 

2 



'* cos 8 <p+ - sin 3 <p + (cos ^ - - sin <p J (l/r)%m ip 
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To show that any value of I > reduces P below its value when 1=0, 
let 



a 


= 


cos 3 tp + - sin 3 ¥> 
P 






( 


1 . x 


b 


= 


I COS <p 


— - sin <p 
P 


c 


= 


MW 










1 - c 
a + 6c 



j sin 2 p 



It suffices to show that 



where the question is in doubt because b may take on negative values. If 
the inequality is to be valid, it is necessary only that (b + a) > 0, that is, 
cos (p > 0. Hence, for the TE modes, only / = needs be considered. For 
this case, the expression for P simplifies to 



(4') 



n > (9) 



n = 0. (10) 







P = 


r 
( 


l 




3 .1.3 

:os <p + - sin <p 
P 


For the 


TM modes, 


there 


is similarly obtained 






P = 


r 

2i 


1 




cos <p -j- - sm <p 
P 






P = 


r 

2^ 


cos <p 




1 . 



cos <p + — sin tp 
2p 



It is easy to show, since cos <p < 1 and sin <p < 1, that both (9) and (10) 
are less than (4'). 

Hence we have shown that, under comparable conditions, i.e., r and p 
constant, the TE Omn modes have higher values of P than any others. 
There is one flaw in the argument, viz., r takes on discrete values and cannot 
be made the same for all modes. It is conceivable, therefore, that for some 
specific values of P, a mode other than the TE Omn can be found which 
gives a smaller W than either of the two "adjacent" TE Omn modes, one 
having a value of r higher, the other lower, than the supposed high-P 
mode. This situation requires further refinement, and hence complication, 
in the analysis; we pass over this point. 

Having so far indicated that the TE Omn modes are the best, our next 
objective is find the best value of m, if possible. 
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By use of the parametric equations (4) and (5), Fig. 2 has been plotted 
for r = 3.83 (TE 01ft modes) and values of n from 1 to 9. This drawing 
shows that, for each discrete value of r, minimum W/P is given by points 
on the "envelope" of the family of curves. 

The standard method of obtaining the envelope is to express W as a 

function of P with n as parameter (r is assumed fixed, for the moment), 

dF ... 

i.e., W = F(P, n), and then set — — 0. However, in this case it is easier 

to express W = G(P, <p) and <p = H(n), whence 

dF &G 6<p 

dn dip dn 

and the envelope is obtained by setting — = provided — ?* 0. We 

proceed, therefore, as follows. 

Assume p is continuous, and solve (4) for p, obtaining: 

-.8 



sin" <p 



r 3 

- — - — cos <p 
2irP 

Now substitute (11) m (5). This gives W as a function of P and ?: 

3 r sin 2 ^ 



W = — 

4tt 2 



2 

COS 



^fe- cosv ). 



(11) 



(12) 



•JTT7 

To solve — = 0, we differentiate and simplify. This yields 
dip 

5 cos 3 ip — 3 cos <p = — p . (13) 

Substituting (13) back into (11) yields 

p = y^ (to 

3 cos 3 ip 

The situation so far is that, with P and r assigned, W lies on the en- 
velope and is a minimum when <p satisfies (13); p is then given by (14). 
Obviously, for (13) to hold, it is necessary that t 

To obtain the best value of r, the procedure is to differentiate W min with 
respect to r, assuming now that r is continuous, and examine for a mini- 
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mum. We can, however, first differentiate (12) by setting 

dw = dw , aw fy 

dr dr dip dr 

dW 
and then substitute from (13). However, when (13) is satisfied, — = 0. 

dip 

This process yields 

dW = r 2 (2 - 3 cos 2 ip) 
dr x 2 9 sin 2 <p cos 6 ip 

dW dW 

This shows — to be positive, when cos 2 *? < f . Hence — = corresponds 

to a maximum, rather than a minimum.* If cos 2 ip < f, that is, ip > 35°16', 
then r should be as small as possible. The smallest r is 3.83, for the TE 
01« modes. For r = 3.83, and ip > 35°, from (13) there is obtained P > 
0.75. 

s 

The analysis thus indicates that, for values of'P = Q- greater than 0.75, 

A 

the TE Olw mode yields the smallest ratio W/P or V/Q. 

An interesting and simple relation between fa and R for minimum W/P 
can easily be derived from the foregoing equations. Substitute (14) back 
into (6), thereby obtaining 

* = ,-^V- (is) 

3a cos 8 <p 

Now use (7) with (15) to eliminate cos <p, replace k by 27r/X, and r by 3.83, 
its numerical value for the TE 01« modes. This gives 

XT* = 223 

or by substituting X = - , c = 3 X 10 1 , 

(fa) 2 R = 20.1 x 10 20 . 

This useful relation was first discovered by W. A. Edson. 

Some further discussion is of interest. It is realized that a number of 
points have not been taken care of in a manner entirely satisfactory mathe- 
matically, but nevertheless important practical results have been obtained. 
As an example, since p and r can assume only discrete values, there are 

* It is for this reason that the determination of the stationary values of W(r, p, <p), 
subject to the constraint P(r, p, ip) = constant, by La Grange multipliers fails to yield 
the desired least value of W/P. 
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specific situations where some mode other than the TE 01m gives a smaller 
W/P. For example, it may be shown that for P between 0.97 and 1.14 
the TE 021 mode yields a smaller W than the TE 013 or TE 014 modes. 
However, the margin is small, and fo» larger P, the TE 02« modes become 
progressively poorer. 

Limitation on Frequency Range of Tunable Cavity as Set 
by Ambiguity 

In the design of a tunable cylindrical resonant cavity intended for use 
in the TE 01m mode, the requirements on Q may dictate a diameter large 
enough to sustain TE 02m' or TE 03m' modes. Also, tbe range of variation 
of cavity length may be such that the TE 01 (m + 1) mode is supported. As 
the cavity is required to tune over a certain range of frequency, the maximum 
frequency range possible in the TE 01m mode without interference from the 
TE 01 (m + l)t or any TE 02 or TE 03 modes is of interest. The interference 
from the TE 01 (m + 1) limits the useful range of the TE 01m by the presence 
of extraneous responses at more than one dial setting for a given frequency 
or more than one frequency for a given dial setting. In applications so far 
made, it has been possible to eliminate extraneous responses from the TE 02 
and TE 03 modes, but crossings of these modes with the main TE 01m mode 
have not been permitted. No designs have had diameters sufficiently large 
to support TE 04 modes. 

The desired relations are easily obtained by simple algebraic manipula- 
tion of equation (1). For simplicity in presentation of the results, we in- 
troduce some symbols applicable to this section only: 



=w b - 1 a - 2w * i& 



A = 

L * J L^J 

A = value of A for TE 01m modes = 13.371 X 10 20 

/ = A/A 

Xo = (a/L) 2 at low frequency end of useful range of TE 01m mode 

maximum/ 

F = frequency range ratio = — —. -. 

minimum / 

The values of A and / depend upon the interfering mode under considera- 
tion. For the TE 02m modes, A = 44.822 X 10 20 , / = 3.3522. 
The two typical cases of interest are shown on Fig. 3. For case I, ara- 

t It is easy to show that the extraneous response from the TE 01 (w — 1) mode is not 
limiting. The proof depends on the inequality n* > (» + 1) (w — 1). 
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TE 01 



Fig. 3— Mode chart illustrating types of interference with TE 01m mode, 
biguity from TE 01(» + 1) mode, it is found that 



F 2 = 



A + B(n + 1) 2 * 
A +Bn 2 Xo 



Curves of F for this case are shown on Fig. 4. 
The maximum value of F is obtained when x = <» and is 



F max 



n + 1 
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Fig. 4 — Curves showing maximum value of frequency ratio without interference from 
TE 01 (« + 1) mode (case I of Fig. 3). 



Table I. — Case II: Maximum Frequency Range Ratio, F, for TE Oln Mode when Limited 
by Mode Crossings with TE 02m and TE 02(m+l) Modes. 





» = 3 


n = 4 


n = 


12 




F 


(<^>min 


F 


("/£)min 


F 


( a / £ >min 


1 


1.198 


1.323 


1.086 


0.966 


1.008 


0.313 


2 






1.242 


1.080 


1.013 


0.316 


3 










1.019 


0.322 


4 










1.027 


0.331 


5 










1.037 


0.343 


6 










1.051 


0.360 


7 










1.071 


0.384 ■ 


8 










1.104 


0.418 


9 










1.168 


0.471 


10 










1.345 


0.564 
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For case II, range limited by mode crossings, it is found that 
A - A 



x = 
F 2 = 



B(n* - «' 2 ) 

(n 2 - n' 2 )[n 2 t - (n' + l) 2 ] 
(w 2 / - «' 2 )[w 2 -(«' + l) 2 ]' 



Some values for this case are given in Table I. 

The formulas above are general and may be used for any pair of mode 
types by using the appropriate values for A and /. 

Tire Elliptic Cylinder 

In the design of high Q circular cylinder cavity resonators operating in 
the TE Oln mode, it is desirable to know how much ellipticity is tolerable, 
so that suitable manufacturing limits may be set. The elliptical wave 
guide has already been studied, notably by Brillouin 12 and Chu, ia but the 
results are not in suitable form or of adequate precision for the present 
purposes. More recently tables 14 have become available which permit the 
calculation of some of the properties of the elliptical cylindrical resonator. 

The elliptical cavity involves Mathieu functions, which are considerably 
more complicated than Bessel functions. 16 The tables give the numerical 
coefficients of series expansions, in terms of sines, cosines, and Bessel func- 
tions, of the Mathieu functions up to the fourth order. These tables have 
been used for the calculation of some quantities of interest in connection 
with elliptical deformations of a circular cylinder in the TE Oln mode. 

The Ellipse 

All mathematical treatments of the ellipse (including the tables men- 
tioned above) use the eccentricity, e, as the quantity describing the amount 
of departure from the circular form. The eccentricity is the ratio 

distance between foci 

e = : : . 

major axis 

This is not a quantity subject to direct measurement, hence we here in- 
troduce and use throughout the ellipticity, E, defined as 

_ difference between major and minor diameters 
major diameter 

It is clear that the ellipticity is easily obtained directly. 

Again, many results are given in terms of the major diameter. Since we 
are interested in deformations from circular, and in such deformations the 
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perimeter remains constant, while the major diameter changes, we have 
expressed our results in terms of an average diameter, defined as 

_ perimeter 

7T 

Figure 5 shows the ellipse and various relations of interest. 

Y 



P=PERIMETER 




2 3 (MAJOR AXIS) 

e=ECCENTRICITY=-^ 

E=ELLIPTICITY=^|p 

A= AREA = TTab 

D = "AVERAGE" D!AMETER = ;£ 



b=aVi-e2 = a (,_ E ) 
A=7ra 2 ~Vi-e 2 =7ra2 (i-e) 



Fig. 5— The ellipse 

Elliptic Coordinates and Functions 

The elliptic coordinate system is shown on Fig. 6. Following Stratton, 15 
we have used £ in place of the table's z, since we wish to use z as the coor- 
dinate along the longitudinal axis. Stratton also uses tj = cos tp as the angu- 
lar coordinate; this is frequently convenient. 

Analogous to cos (6 and sin (6 in the circular case, there are even and 
odd* angular functions, denoted by 

e S/(c, cos <p) and "Sf(c, cos <p) 



which reduce to cos Id and sin 16 respectively when c 
are even and odd* radial functions, denoted by 

'J t (c, h) and V/(c, £) 

* For t = 0, only even functions exist. 



0. Similarly, there 
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which both reduce to Jf(kip) when c — ■> 0. In the above, c is a parameter 
related to the ellipticity.* The tables do not give values of the functions, 
but rather give numerical coefficients 

Di and F{ 

of expansions in series of cosine, sine and Bessel functions, which permit one 
to calculate the elliptic cylinder functions. The coefficients, of course, 



60 **/, 




Fig. 6 — Elliptic coordinate system 



depend on the parameter c\ the largest value of c in the tables is 4.5, which 
corresponds to an ellipticity of 39% in a cylinder operating in the TE 01/? 
mode.** For this case, Bessel functions up to Jn(x) and Jn(x) are needed 
for calculating the radial function. It is clear that calculations on elliptic 
cylinders have not been put on a simple basis. 

* Not to be confused with c = velocity of electromagnetic waves; the symbol c is 
here carried over from the published tables. 

** An ellipticity of 39% means that the difference between maximum and minimum 
diameters is 39% of the maximum diameter. For a given c, the ellipticity depends on 
the mode- 
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Field Equations 

The equations for the fields are easily obtained from section 6.12 of 
Stratton's book, and are given in Table II, which is self-explanatory, except 
for the quantity c, which we now proceed to discuss. 

Resonant Frequencies 

The elliptic cylinder has the major diameter, 2a, and the focal distance, 

a 
2c . The equation of its surface is then expressed by £ = - = a. On 

this surface, E„ must vanish. This requires that '°J/,(c, a) = for TE 
modes and that '°J/,(c, a) = for TM modes. The series expansions are 
ui terms of c£ as variable. Let ca = r' lm or r £m be the roots of the above 

equations. Then - = - (dropping the subscripts C, m). Now, in working 
Co a 

out the solution of the differential equations, it turned out that c = c h. 
Here h is one component of the wave number, k]. Hence h = - . Further- 

more, the eccentricity is e = - = - . The indicated procedure is: 1) choose 

' J a r 

a value of c; 2) find the various values of r for which the radial function or 
its derivative is zero; 3) then calculate the corresponding eccentricity and 
resonant frequency. Notice that for a given value of c, the values of r 
will depend on the mode, and hence so will the eccentricity. 

We now wish to express our results in terms of the ellipticity and the 
average diameter. To convert eccentricity to ellipticity, we use 

E = 1 - Vl - e 2 - 

The perimeter of the ellipse is given by P = AaE(e) where E{e) is the com- 
plete elliptic integral of the second kind. ft 
In terms of the average diameter we find 

_ 2|"2r£(e)1 

or calling the quantity in brackets s, h = j: ■ This is now in the same form 
as ki for a circular cylinder of diameter D. The quantity s is the recipro- 
cal of Chu's -^. 
t It is recalled that 

k = — = Vk* + k\ ; *. = - ; h = — . 

A » J a L 

ft This is tabulated as E(a) in Jahnke & Emde, p. 85, with a = sin^e. 
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We have calculated and give in Table III values of r, e, E and s for several 
values of c and for a few modes of special interest. For three cases, e TE 01, 
"TM 11 and "TM 11, we have determined an empirical formula to fit the 
calculated values of s. These are also given in Table III. 



TE Modes 



TABLE II. Elliptic Cylinder Fields 



E * = ~ k A/- ~ ^i v)Jf(c, & sin k» > cos ut 

E„ = k i/f! ^i^- 1 Sg(c, ii)J'((c, Q sin k,z cos ut 
y/ti — 1 

Ht = h S((c, rj)J[(c, Q cos faz sin ut 

1 

77. = £3 S((c, tj)7/(c, £) cos £ 3 z sin ut 

1 

n z = k'fSf(c, «j)/^(c, *) sin k a z sin ut 



TM Modes 



\A 2 — 1 

E$ = —£3 S((c, rj)J((c, £) sin £ 3 2 cos ut 

1 

E„ = —k 3 S((c, i))J/(c, £) sin k 3 z cos ut 

9 

E z = k\ S((c, r))J((c, cos k 3 s cos ut 

n£= -k A /- " S£(c, i|)//(e, f) cos A, 3 sin ut 

/" \/fc2 _ J 

H v = k A/ - — •fyCc, r))J((c, £) cos fej2 sin ut 

V H q 

Notes: 
Derivatives are with respect to £ and tj. 
5/ and J / carry prefixed superscripts, c or 0, since they may be either even or odd. 

q = c Q V$ 2 — »? a c = coki 

ky = T -^ k, = "y V = k'[ + k] 

a L 

2co is distance between foci of ellipse. 
a is the semi major diameter of the ellipse. 
r s „, is the value of c$ that makes 

Jf(cji) = for TM modes 
Jfc£) = for TE modes. 
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TABLE III Root Values of Radial Elliptic Cylinder Functions 



Mode 


c 


r 


e 


E 


s 


'TE 01 





3.8317 








3.8317 




0.2 


3.8343 


0.05216 


0.001361 


3.8317 




0.4 


3.8423 


0.10410 


0.005434 


3.8318 




0.6 


3.8558 


0.15561 


0.012181 


3.8324 




0.8 


3.8753 


0.20643 


0.021539 


3.8337 




1.0 


3.9015 


0.25631 


0.033406 


3.8366 




1.2 


3.9349 


0.30496 


0.047636 


3.8417 




1.4 


3.9763 


0.35209 


0.064033 


3.8500 




1.6 


4.0264 


0.39738 


0.082346 


3.8624 




2.0 


4.154 


0.4814 


0.12351 


3.902 




3.0 


4.634 


0.6474 


0.2378 


4.101 




4.0 


5.29 


0.756 


0.346 


4.42 




4.5 


5.66 


0.795 


0.393 


4.62 





.v = 


3.8317 + 4.33 E 2 + 1.9£ 3 






'TM 11 





3.8317 








3.8317 




0.2 


3.8330 


0.05218 


0.001362 


3.8304 




0.4 


3.8370 


0.10425 


0.005449 


3.8265 




0.6 


3.8436 


0.15610 


0.012259 


3.8201 




0.8 


3.8532 


0.20762 


0.021791 


3.8113 




1.0 


3.8658 


0.25868 


0.034036 


3.8003 




1.2 


3.8818 


0.30913 


0.048981 


3.7874 




1.4 


3.9015 


0.35884 


0.066599 


3.7727 




1.6 


3.9253 


0.40761 


.0.086844 


3.7568 




4.5 


5.13 


0.878 


0.520 


3.91 



j = 3.8317 - 0.96E + LIE* 



°TM 11 





3.8317 








3.8317 




0.2 


3.8356 


0.05214 


0.001361 


3.8330 




0.4 


3.8474 


0.10397 


0.005419 


3.8370 




0.6 


3.8670 


0.15516 


0.012111 


3.8436 




0.8 


3.8944 


0.20542 


0.021326 


3.8530 




1.0 


3.9298 


0.25446 


0.032918 


3.8654 




1.2 


3.9731 


0.30203 


0.046701 


3.8809 




1.4 


4.0243 


0.34788 


0.062462 


3.8997 





s = 


= 3.8317 + 0.95£ + 2.2E 2 






'TE 22 





6.706 








6.706 




0.4 


6.712 


0.0596 


0.00178 


6.706 




0.8 


6.729 


0.1189 


0.00709 


6.705 




1.2 


6.756 


0.1776 


0.01590 


6.702 




1.6 


6.788 


0.2357 


0.02817 


6.693 




2.0 


6.826 


0.2930 


0.04389 


6.677 


°TE 22 





6.706 








6.706 




0.4 


6.712 


0.0596 


0.00178 


6.706 




0.8 


6.730 


0.1189 


0.00709 


6.706 




1.2 


6.762 


0.1775 


0.01587 


6.708 




1.6 


6.810 


0.2350 


0.02799 


6.715 




2.0 


6.877 


0.2908 


0.04323 


6.729 
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Mode 


c 


r 


e 


E 


s 


•TE 32 





8.015 








8.015 




0.4 


8.020 


0.0499 


0.00124 


8.015 




0.8 


8.035 


0.0996 


0.00497 


8.015 




1.2 


8.059 


0.1489 


0.01115 


8.014 




1.6 


8.093 


0.1977 


0.01974 


8.013 




2.0 


8.135 


0.2459 


0.03070 


8.010 


•TE 32 





8.015 








8.015 




0.4 


8.020 


0.0499 


0.00124 


8.015 




0.8 


8.035 


0.0996 


0.00497 


8.015 




1.2 


8.060 


0.1489 


0.01115 


8.015 




1.6 


8.097 


0.1976 


0.01972 


8.018 




2.0 


8.146 


0.2455 


0.03061 


8.022 


•rjf oi 





2.4048 











0.2 


2.4090 


0.08302 








0.4 


2.4216 


0.16518 








0.6 


2.4431 


0.24559 








0.8 


2.4739 


0.32337 








1.0 


2.5149 


0.39762 






TEll 





1.8412 











0.2 


1.8416 


0.10860 








0.4 


1.8430 


0.21704 








0.6 


1.8452 


0.32516 








0.8 


1.8484 


0.43280 








1.0 


1.8527 


0.53975 







Notes: 

Superscripts e and o on mode designation signify even and odd. 

c is parameter used in the Tables (Stratton, Morse, Chu, Hutner, "Elliptic Cylinder 

and Spheroidal Wave Functions") 
r is the value of the argument which, for TM modes, makes the radial function zero 

and, for TE modes, makes its derivative zero. 
e is the eccentricity of the ellipse; 

_ distance between foc i 

major diameter 

E is the ellipticity of the ellipse; 

_ difference between major and minor diam. 

major diameter 

5 is the root value, referred to the "average diameter"; it is related to r by: 



r 
s = - 



perimeter 



it major diameter 

The quantity j is also related to the cutoff wavelength in an elliptical wave guide 
according to: 

perimeter of guide 



cutoff wavelength 



Resonator Q 



Although the calculation of the root values is straightforward and not 
overly laborious, the same cannot be said for the integrations involved in 
the determination of resonator Q. The procedure is obvious: The field 
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equations are given; it is only necessary to integrate EHr over the volume 
and H 2 da over the surface and get Q from 

2 / R2dT 
Q = \ L (16) 

S J B>d„ 

with 5 = skin depth, a known constant. Unfortunately the integrations 
cannot at present be expressed in closed form. A numerical solution can 
be obtained by a combinaticn of integration in series and of numerical 
integration. 

The calculations have been made for the 'TE 01 mode with c = 2.0, for 
which r = 4.154. This value of c corresponds in this case to an ellipticity 
of about 12%; in a 4" cylinder this would amount to 1/2" difference between 
largest and smallest diameters. Evaluation* of the integrals yields: 



H 2 dr = 12.307 k\L + 12.294 k\L 
v 



H 2 d<r = 49.228 k\ + 0.1619 kyk\L + 6.6847 k\L 

s 

7.804 , , irn , . . £ „ 

Substituting k\ = —=r- and k 3 = — , one obtains, finally 

/ 1 + 0.1622 n 2 R 2 \ 

Q8 - 0.471 D ^ + Q Q039 n2R2 + Q 1529 n2R3 J • 



For a circular cylinder, 



Q c 8 = 0.5 D 



2tj2 



1 + 0.1681 h'R 
f+ 0.1681 n*R 3 



Comparison of these two formulas for Qd shows that the losses in the end 
plates (nrR 3 term) are less with respect to the side wall losses in the ellip- 
tical cylinder. The net loss in Q8, as described by the reduction in the mul- 
tiplier frcm 0.5 to 0.471, is thus presumably ascribable to an increase in side 
wall losses (stored energy assumed held constant). The additional term 
in n 2 R- in the denominator is responsible for the difference in the attenuation- 
frequency behavior of elliptical vs circular wave guide as shown by Chu, 
Fig. 4. Incidentally, these results agree numerically with those of Chu. 

* Numerical integration was by Weddle's rule; intervals of 5° in <p and 0.1 in x were used. 
The calculations were made by Miss F. C. Larkey. 
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Corresponding expressions for the resonant wavelength are 
ttD 0.805 D 



X = 



X e = 



0.820 D 
VI + 0.1681 « 2 U 2 " 



Vl + 0.1622 w 2 2? 2 



As an example, take n = 1, R = 1, then 

(Circular) &S = 0.500 D X c = 0.759 D 

(Elliptical) Q8 = 0.473 D X = 0.747 D 

Ratio = 0.946 Ratio = 0.984. 

Conclusions 

The mathematics of the elliptic cylinder have not yet been developed to the 
point where the design of cavities of large ellipticity could be undertaken. 
On the other hand, sufficient results have been obtained to indicate that the 
ellipticity in a cavity intended to be circular, resulting from any reasonable 
manufacturing deviations, would not have a noticeable effect on the reso- 
nant frequencies or Q values, at least away from mode crossings. 

Full Cylindrical Coaxial Resonator 

The full coaxial resonator has been of some interest because of various 
suggestions for the use of a central rod for moving the tuning piston in a 
TE 01m cavity. 

The cylindrical coaxial resonator, with the central conductor extending 
the full length of the resonator, has modes similar to the cylinder. In 
fact, the cylinder may be considered as a special case of the coaxial. The 
indices /, m, n have much the same meaning and the resonant frequencies 
are determined by the same equation (1). However, now the value of r 
depends in addition (see Fig. 1) upon 77, where 

_ diameter inner conductor _ b 
diameter outer conductor a 

The problem now arises of how best to represent the relations between 
/, o, b and L. The r's depend on 77; so one possibility is to determine their 
values for a given 17 and then construct a series of mode charts, one for each 
value of 77. 

A more flexible arrangement is to plot the values of r vs 77 and allow 
the user to construct graphs suitable for the particular purpose in hand. 
An equivalent scheme has been used by Borgnis. 18 

It turns out that as r\ — > 1, r(l — 77) — » mir, for the TM modes and the 
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TE Oro« modes, and r(l — ij) — » (m — 1)tt for all other TE modes. For 
the former modes, r becomes very large as r; — > 1, that is, as the inner con- 
ductor fills the cavity more and more, the frequency gets higher and higher. 
For the TE l\n modes, however, as the inner conductor grows, the fre- 
quency falls to a limiting value. This is discussed in more detail by 
Borgnis. 18 

Figure 7 shows r(l — tj) vs -q, for a few of the lower modes; the scale for r\ 
between 0.5 and 1.0 is collapsed since this region does not appear to be of 
great engineering interest. A different procedure is used for the roots of 
the TE fin modes. Figure 8 is a direct plot of r vs 77 for a few of the lower 
modes. In this case, r — > I as ij — > 1 . 

Distribution of Normal Modes 

The calculation of the distribution of the resonant modes for the coaxial 
case follows along the lines of that for the cylinder, as given previously. 
The difference lies in the distribution of the roots r, which now depend upon 
the parameter r,. The determination of this latter distribution offers 
difficulties. There is some evidence, however, that the normal modes will 
follow, at least to a first approximation, the same law as the cylinder, viz.: 

V 

N = 4.4 73 

Ao 

with some doubt regarding the value of the coefficient. 
Q - in Coaxial Resonator 

A 

The integrations needed to obtain this factor are relatively straightfor- 
ward, but a little complicated. The final results are given in Fig. 1. 

The defining equation is (16); the components of H are givsn in Fig. 1. 
The integrations can be done with the aid of integrals given by McLachlan 17 
and the following indefinite integral : 

/[/ 2 i*W + #(*)]*& 

which can be verified by differentiation, remembering that y = Z((x) is a 
solution of y" + - y' + ( 1 A y = 0. 
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An investigation needs to be made of the behavior of the formulas as 
77 — * before any conclusion may be drawn regarding their blending 
into those for the cylinder. For TE modes with I = 0, the term involving 

— disappears, hence no question arises. Consider then I > 0, and let 

V 

x = tjr for the discussion following. From expansions given in McLachlan, 

it is easy to show that, for small * 
(( — 1)' /2V 



Since, from Fig. 1, 



■K \x) X 






X 1 

'*> - 2T7: 


f((x) « - 


*r* 


\*{t - 1) 1 



#(') F/(ifr) Fftz) 
it is found, upon substitution of the approximations given above: 

2x* 

That is, Z((x) ~ a^ and hence — > as x — > 0. Furthermore Z;(r) remains 
finite as tj -► 0. Hence H ~ X H and - ~ a 2 * -1 . Therefore, for / > 0, 

- -* as rj -> 0. 

Hence, the expression for Q - for the coaxial structure reduces to that for 

A 

the cylinder, for any value of I, in the TE modes. 
For the TM modes, and for C > 0, an entirely similar argument shows 

that H' remains finite as ij — > 0. Hence, the expression for Q - for these 

A 

modes also reduces to that for the cylinder. 

For the TM modes, and with / = 0, we have 

Z' (x) - -/,(*) + /.(») Dg . 

For x —> 0, 7](.r) — > and Jo(x) — * 1, hence for small x, 

Zo(x) ~ m ■ 
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Now substitute the approximate values of the Y for small x. The result is 

1 



Zo(-r) r 
Since Z' (r) is finite, it follows that 



***\ 



* ( l0g t) 



,217' 



and it is easily shown that rjH ' — > oo as r\ -* 0. On the other hand, 77 2 # 
as j; — > 0. Hence, Q - -> as ?? ->• 0. On the other hand, for ?j = 0, a 

A 
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Fig. 9— Coaxial resonator. TE Oil mode Contour lines of Gr- 



it 
perfect cylinder exists whose Q - is not zero. It is concluded that the ex- 
pression for Q - does not apply for small r? for the TM modes with t = 0. 

A 

Thus it is seen that the expressions for the factor (Q -) reduce to those 

A 

given for the cylinder, when -q = 0, except for TM modes with C = 0. 

For these latter cases, the factor approaches zero as t/ approaches zero, 

because rjH' increases without limit, This means that an assumption 
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Fig. 10 — Coaxial resonator. TE 111 mode Contour lines of Q- 
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which was made in the derivation of the Q values is not valid for small rj; 
that is, the fields for the dissipative case are not the same as those derived 
on the basis of perfectly conducting walls. 

The expressions for the factor are rather complicated, as it depends on 
several parameters. When a given mode is chosen, the number of param- 
eters reduces to two, 17 and R. Contour diagrams of Q - vs rj and R are 

A 

given on Figs. 9, 10, 11 and 12 for the TE Oil, TE 111, TM Oil and TM 111 
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Fig. 12— Coaxial resonator. TM 111 mode Contour lines of Q- 



modes. As mentioned above, the true behavior of Q - for the TM Oil 

mode for small 77 is not given by the above formula, so this contour diagram 
has been left incomplete. 

Fins in a Cavity Resonator 

The suppression of extraneous modes is always an important problem 
in cavity design. Among the many ideas advanced along these lines is the 
use of structures internal to the cavity. 

It is well known that if a thin metallic fin or septum is introduced into a 
cavity resonator in a manner such that it is everywhere perpendicular to 
the E-lines of one of the normal modes, then the field configuration and 
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frequency of that particular mode are undisturbed. For example, Fig. 13 
shows the £-lines in a TE 11» mode in a circular cylinder. If the upper 
half of the cylinder wall is replaced by a new surface, shown dotted, the 
field and frequency in the resulting flattened cylinder will be the same as 



NEW SURFACE PERPENDICULAR 
TO E- LINES LEAVES REST OF 
FIELD AND FREQUENCY UNALTERED 




Fig. 13 — E Lines in TE 1 In mode 



^ORIGINAL CYLINDER 



i-' 




Fig. 14 — "TE Oln" mode in half-cylinder 



before. Indeed, they will also be the same in the crescent-shaped resonator 
indicated in the figure. 

Except for isolated cases, all the other modes of the original cylinder will 
be perturbed in frequency since the old fields fail to satisfy the boundary 
conditions over the new surface. Furthermore, if the original cylinder was 
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circular, its inherent double degeneracy will be lost and each of the original 
modes (with minor exceptions) will split into two. 

Although the frequency and fields of the undisturbed mode are the 
same, the Q is not necessarily so. For example, Fig. 14 shows a "TE 01« 
mode" in a half cylinder.* 



It is easy to calculate Q - for this case. The result is 

A 

n 8 __r_ (1 + P 2 RV 12 

V \ 2t 1 + p 2 R* + Kx + KiP*R* 



(17) 



in which 



K x = 1.290 K 2 = 0.653. 



Here K\ and K% are constants which account for the resistance losses in 
the flat side. For the full cavity, shown dotted in Fig. 14, eq. (17) holds 
with K\ = K 2 = 0. If the circular cavity has a partition extending from 
the center to the rim along the full length, (17) holds with the values 
of K\ and Ki halved. If a fin projects from the rim partway into the in- 
terior, still other values of K x and Ki are required. It is a simple matter 
to compute these for various immersions; Fig. 15 shows curves of K\ and 
K 2 . The following table gives an idea of the magnitudes involved: 

mode: TE 0,1,12 R = 0.4 



Fin,%o 


4 


Ratio 


0% 


2.573 


1.0 


10 


2.536 


.985 


20 


2.479 


.965 


so 


2.04 


.79 


100 


1.47 


.57 



The question now is asked, "Suppose a longitudinal fin were used, small 
enough to cause only a tolerable reduction in the Q. Would such a fin 
ameliorate the design difficulties due to extraneous modes?" 

Some of the effects seem predictable. All modes with t > will be split 
to some extent, into two modes of different frequencies. Consider the 
TE Yin mode, for example. There will be one mode, of the same frequency 
as the original whose orientation must be such that its E-lines are perpendicu- 
lar to the fin. The Q of this mode would be essentially unchanged. There 
will be a second mode, oriented generally 90° from the first, whose £-lines 
will be badly distorted (and the frequency thereby lowered) in the vicinity 

* Solutions for a cylinder of this cross-section are known and all the resonant fre- 
quencies and Q values could be computed, if they had any application. 
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of the fin. It would be reasonable to expect the Q of this mode to be appre- 
ciably lowered because of the concentrated field there. If two fins at 90° 
were present, there would be no orientation of the original TE Yin mode 
which would satisfy the boundary conditions. In this case both new modes 
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Fig. 15 — Constants for calculation of Q of TE Olw mode in cylinder with longitudinal 



would be perturbed in frequency from the original value. If both fins were 
identical, the perturbations would be equal and a double degeneracy ensue. 
Similar effects would happen to the other types of modes. 

The major advantage derivable from such effects would appear to be in 
extraneous transmissions. The fin serves to orient positively the fields in 



442 BELL SYSTEM TECHNICAL JOURNAL 

the cavity, and the input and output coupling locations can then be appro- 
priately chosen. On the basis that internal couplings are responsible for 
mode crossing difficulties, one might hazard a guess that a real fin would 
increase such couplings. 

Another application of fins might be in a wave guide feed in which it is 
desired to establish only a TE Om wave. In this case, Q is not so important 
and larger fins can be used. If these extended virtually to the center and x 
of them were present (with uniform angular spacing) all types of wave trans- 
mission having I less than x/2, x even or ( less than x, x odd, would be sup- 
pressed. This use of fins is an extension of the wires that have been 
proposed in the past. 

Conclusion 

It is hoped that the foregoing, which covers some of the theoretical work 
done by the author during the war, will be of value to other workers in 
cavity resonators. There is much that needs to be done and hardly time 
for duplication of effort. 
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